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Abstract. In this work we introduce a new family of 14-steps linear multistep 
methods for the integration of the Schrodinger equation. The new methods are 
phase fitted but they are designed in order to improve the frequency tolerance. This 
is achieved by eliminating the first derivatives of the phase lag function at the fitted 
frequency forcing the phase lag function to be 'flat' enough in the neighbor of the 
fitted frequency. The efficiency of the new family of methods is proved via error 
analysis and numerical applications. 
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1. Introduction 



The numerical integration of systems of ordinary differential equations 
with oscillatory solutions has been a subject of research during the past 
decades. This type of ODEs is often met in real problems encounter in 
computational chemistry, like the Schrodinger equation. For problems 
having highly oscillatory solutions standard methods with unspecial- 
ized use can require a huge number of steps to track the oscillations. 
One way to obtain a more efficient integration process is to construct 
numerical methods with an increased algebraic order, although the 
implementation of high algebraic order meets several difficulties like 
resonances [3] . 
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On the other hand, there are some special techniques for optimizing 
numerical methods. Trigonometrical fitting and phase-fitting are some 
of them, producing methods with variable coefficients, which depend on 
V = Loh, where lo is the dominant frequency of the problem and h is the 
step length of integration. More precisely, the coefficients of a general 
linear method arc found from the requirement that it integrates exactly 
powers up to degree p + 1. For problems having oscillatory solutions, 
more efficient methods are obtained when they are exact for every linear 
combination of functions from the reference set 

{l,x,...,x^e±'^^...,x^e±'^-} (1) 

This technique is known as exponential (or trigonometric if = icu) 
fitting and has a long history [2], [1]. The set (1) is characterized by 
two integer parameters, K and P . The set in which there is no classical 
component is identified hy K = — 1 while the set in which there is no ex- 
ponential fitting component (the classical case) is identified by P = — 1. 
Parameter P will be called the level of tuning. An important property 
of exponential fitted algorithms is that they tend to the corresponding 
classical ones when the involved frequencies tend to zero, a fact which 
allows to say that exponential fitting represents a natural extension of 
the classical polynomial fitting. The examination of the convergence 
of exponential fitted multistep methods is included in Lyche's theory 
[1] . There is a large number of significant methods presented with high 
practical importance that have been presented in the bibliography. The 
general theory is presented in detail in [5]. 

Considering the accuracy of a method, when solving oscillatory prob- 
lems, it is more appropriate to work with the phase-lag, rather than 
the principal local truncation error. Wc mention the pioneering paper 
of Brusa and Nigro [4], in which the phase-lag property was intro- 
duced. This is actually another type of a truncation error, i.e. the angle 
between the analytical solution and the numerical solution. On the 
other hand, exponential fitting is accurate only when a good estimate 
of the dominant frequency of the solution is known in advance. This 
means that in practice, if a small change in the dominant frequency is 
introduced, the efficiency of the method can be dramatically altered. 
It is well known that for equations similar to the harmonic oscillator 
the most efficient exponential fitted methods are those with the highest 
tuning level. A lot of significant work has been made during the last 
years in this field, mainly focusing for obvious reasons in the solution 
of the Schrodinger equation (see for example [10]- [132]). 

In this paper we present a new family of methods based on the 14- 
step linear multistep method of Quinlan and Tremainc [6]. The new 
methods are constructed by vanishing the phase-lag function and its 
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first derivatives at a predefined frequency. Error analysis and numeri- 
cal experiments show that the new methods exhibit improved charac- 
teristics concerning the solution of the time-independent Schrodinger 
equation. The paper is organized as follows: In section 2, the general 
theory of the new methodology is presented. In section 3, the new 
methods are described in detail. In section 4 the stability properties of 
the new methods are investigated. Section 5 presents the results from 
the numerical experiments and finally, conclusions are drawn in section 
6. 



2. Phase-lag analysis of symmetric multistep methods 

Consider the differential equations 

^ = fit,y), y{to) = yo, y'(to) = y'o (2) 

and the linear multistep methods 

J J 

y^+j = ^^Y^ ^0 fn+j (3) 

j=0 j=0 

where yn+j = y{to + {n + j)h), fn+j = f{to + {n + j)h,y{to + {n + j)h)) 
and h is the step size of the method. With the method (3) , we associate 
the following functional 

J J 
Lih, a, b, yit)) = + j-h)-h^Yl hy\^ + J ' h) (4) 

j=o j=o 

where a, b are the vectors of coefficients aj and bj respectively, and y{t) 
is an arbitrary function. The algebraic order of the method (3) is p, if 

L{h, a, b, y{t)) = Cp+2hP+^y^+^\t) + 0{hP+^) (5) 
The coefficients Cq are given 

Co = Y.j=o O'j 

Ci = E/=oi-«.- _ (6) 

The principal local truncation error (PLTE) is the leading term of (5) 

PLTE = Cp+2/i^+^y(^+^) (t) (7) 
The following assumptions will be considered in the rest of the paper: 
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1. aj = 1, since we can always divide the coefficients of (3) witfi aj. 

2. |ao| + |6o| 7^ 0, since otherwise we can assume that J = J — 1. 

3. J2j=o\^j\ / 0' since otherwise the solution of (3) would be inde- 
pendent of (2). 

4. The method (3) is at least of order one. 

5. The method (3) is zero stable, which means that the roots of the 
polynomial 

J 

p{z) = J2ajz^ (8) 

j=0 

all lie in the unit disc, and those that lie on the unit circle have 
multiplicity one. 

6. The method (3) is symmetric, which means that 

% = bj = bj-j, j = 0(1)J (9) 

It is easily proved that both the order of the method and the step 
number J are even numbers [7]. 

Consider now the test problems 

y"{t) = -cj^yit) (10) 

where a; is a constant. The numerical solution of (10) by applying 
method (3) is described by the difference equation 

J/2 

^ Aj{s'^){yn+j + Vn-j) + Ao{s'^)yn = (11) 

with 

Aj{s^)=ai_. + s^ ■bi_. (12) 
and s = uh. The characteristic equation is then given by 

J/2 

Y,A,{s^){z^ +z-^) + Ao{s') = Q (13) 
i=i 

and the interval of periodicity (0, Sq) is then defined such that for s G 
(0, So) the roots of (13) are of the form 

Zl = e'^^'\ Z2 = e-'^^'\ \zj\ < 1, 3 < j < J (14) 
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where A(s) is a real function of s. The phase-lag PL of the method (3) 
is then defined 

PL = s- X{s) (15) 

and is of order q if 

PL = c- + 0(s«+=^) (16) 

In general, the coefficients of the method (3) depend on some parameter 
V, thus the coefficients Aj are functions of both and v. The following 
theorem was proved by Simos and Williams [8]: For the symmetric 
method (10) the phase-lag is given 

2 e/£i Ms\v)-cos{j-s)+ Ao{s\v) 

We are now in position to describe the new methodology. In order 
to efficiently integrate the Schrodinger equation, it is a good practice 
to calculate the coefficients of the numerical method by forcing the 
phase lag to be zero at a specific frequency. But, since the appropriate 
frequency is problem dependent and in general is not always known, 
wc may assume that we have an error in the frequency estimation. It 
would be of great importance to force the phase-lag to be insensitive to 
this error. Thus, beyond the vanishing of the phase-lag, we also force 
its first derivatives to be zero. 



3. Construction of the new methods 
3.1. Classical Method 

The family of new methods is based on the 14-step linear multistep 
method of Quinlan and Tremaine [6] which is of the form (3 with 
coefficients 



ao = 1 ao = —2 a2 = 2 03 = —1 04 = 05 = uq = ay = 



60 = 

_ 433489274083 28417333297 

^ ~ 237758976000 ^ ~ ~ 4953312000 ^ 

176930551859 , 7854755921 , 

64 = 65 = 06 = 

o^^n^o^^^^ 65228800 



2971987200 
_ 577045151693 

^ ~ 2830464000 



_ 930518896733 

" 39626496000 

146031020287 

825552000 



(18) 
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The PLTE of the method is given by 

152802083671y(^6)/,i6 / x 

PLTE = + 0(h^^] (19) 

2853107712000 V 7 ^ ^ 

3.2. New Methods using Phase Fitting 

The methods that are constructed are named as PF-Di, where: 

— PF-DO: the phase lag function is zero at the frequency v = u) *h. 

— PF-Dl : the phase lag function and its first derivative are zero at 
the frequency v = u * h. 

— PF-D2: the phase lag function and its first and second derivatives 
are zero at the frequency v = u *h. 

— PF-D3: the phase lag function and its first, second and third 
derivatives are zero at the frequency v = lo * h. 



— PF-D4 ■ the phase lag function and its first, second, third and 
fourth derivatives are zero at the frequency v = lo * h. 



— PF-D5: the phase lag function and its first, second, third, fourth 
and fifth derivatives are zero at the frequency v = u *h. 



— PF-D6: the phase lag function and its first, second, third, fourth, 
fifth and sixth derivatives are zero at the frequency v = to * h. 

The coefficients of the methods in the form: 



b\ 



l.nurn 

w ■ 

denum 
"4,nttm 

denum 



bi 



2/num 

y^. — ' 

. denum 
hi 

"5,num 

V' ' 

denum 



3,num 
denum 



denum 



b\ 



^7,num 
hi 

denum 



(20) 



where the coefficients ¥ correspond to the method PF-Di. Since for 
small values of v, the above formulae are subject to heavy cancelations, 
the Taylor expansions of the coefficients have been calculated as 6^. The 
exact formulae of all coefficients are given in appendix. 
The PLTEs of the methods are: 
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PLTE'' 



PLTE' 



PLTE^ = 



152802083671 

2853107712000 
152802083671 

2853107712000 
152802083671 



PLTE-" 



PLTE* 



PLTE"" 



PLTE" = 



2853107712000 

3y(i^)a;2 + y(i6))/ji6 + o(/,i8 

152802083671 



(y^^^^a;^ + 2y^'''>u;' + y^'^))h'*' + 0{h 
(y(^°)a;6 + 3y(i2)^4^ 



2853107712000 
152802083671 



+ 4^(10)^6 + gy(i2)^4 ^ 4^(14)^2 



2853107712000 
+5y(i4)a;2 + y(i6))^i6 + o(/,i8) 



152802083671 



2853107712000 



152802083671 



2853107712000 
+35y(io)a;6 + 21y('^^u^ + 7y(i^)u;2 ^ y(i6))^i6 + ^(^^iS) 



(21) 



4. Stability Analysis 

The stability of the new methods is studied by considering the test 
equation 



(22) 



and the linear multistep method (3) for the numerical solution. In the 
above equation a ^ lo [to \s the frequency at which the phase-lag 
function and its derivatives vanish). By setting s = ah and v = coh 
we get for the characteristic equation of the applied method 



J/2 



(23) 
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where 

Aj{s'',v) = aj_.{v) + ■ bj_.iv) (24) 

The motivation of the above analysis is straightforward: Although the 

coefficients of the method (3) are designed in a way that the phase-lag 
and its first derivatives vanish in the frequency lo, the frequency to itself 
is unknown and only an estimation can be made. Thus, if the correct 
frequency of tlic problem is a we have to check if the method is stable, 
that is if the roots of the characteristic equation lie in the unit disk. For 
this reason we draw in the s — v plane the areas in which the method 
is stable. Figure 1 shows the stability region for the six methods (the 
classical one, the phase fitted one and those with first, second, third, 
fourth, fifth and sixth phase lag derivative elimination). Note here that 
the s-axis corresponds to the real frequency while the v-axis corre- 
sponds to the estimated frequency used to construct the parameters of 
the method. 



5. Numerical Results 

The radial Schrodinger equation is given by: 

y"i^) = (^^^ + ^(^) - E)y(^)) (25) 

where ^-^^^^^ is the centrifugal potential, V{x) is the potential, E is the 

Energy and W{x) = ^'^^^P + V{x) is the effective potential. It is valid 
that linix^ooy {x) = and therefore lirrix^ooW {x) = 0. We consider 
that > and we divide the interval [0, +co) into subintervals [a,, hi\ 
so that W{x) can be considered constant inside each subinterval with 
value Wj. The problem (25) can be expressed now by the equations 

y'l = {W, - E)yi (26) 

whose solution are 

Vi {x) = e^^" + 5^6-^^^) (27) 

with Ai,Bi e R. We will integrate problem 25 with Z = at the interval 
[0, 15] using the well known Woods-Saxon potential: 

Vix) = —— + J7-—T2, q = e a (28) 
1 + q (1-l-g)^ 

where uq = —50, a = 0.6, xq = 7, ui = — ^ and with boundary 
condition y(0) = 0. The potential V{x) decays more quickly than 
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^^^^P , SO for large x (asymptotic region) the Schrodinger equation 
(25) becomes 

y"{x) = (^^^ - E)y{x)^ (29) 

The last equation has two linearly independent solutions kxji (kx) and 
kxni{kx), where ji and n; are the spherical Bessel and Neumann func- 
tions and k = ^ ^(|+^) When x ^ oo the solution takes the 
asymptotic form 

y{x) ~ Akxji(kx) - Bkxni{kx) , . 

D[sin{kx — TT^) + tan{5i)cos{kx — TT^)], 

where 5i is called the scattering phase shift and it is given by the 
following expression: 

' y{xiJ^xC{xi)-y{xi)C{xi+i) 

where S(x) = kxji{kx) and C{x) = kxni{kx) and Xi < Xj+i and both 
belong to the asymptotic region. Given the energy, wc approximate the 
phase shift, the accurate value of which is ^ for the above problem. 
We will use three different values for the energy: i) 989.701916, ii) 
341.495874 and iii) 163.215341. As for the frequency uj we will use the 
suggestion of Ixaru and Rizea [9] : 

_/ ^/E^, [0,6.5] 
""-{VE, xG[6.5,15] 

The results arc shown in figures 2, 3 and 4. It is clear that the 
accuracy increases as the number of the eliminated derivatives of the 
phase lag function increases. 



6. Conclusions 

We have presented a new family of 14-stcps symmetric multistep nu- 
merical methods with improved characteristics concerning the integra- 
tion of the Schrodinger equation. The methods were constructed by 
adopting a new methodology which, except for the phase fitting at a 
predefined frcqiicncy, it eliminates the first derivatives of the phase lag 
function at the same frequency. The result is that the phase lag function 
becomes less sensitive on the frequency near the predefined one. This 
behavior compensates the fact that the exact frequency can only be 
estimated. Experimental results demonstrate this behavior by showing 
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that the accuracy is increased as the number of the derivatives that are 
ehminated is increased. 
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Appendix 



Method PF-DO: 

bi,num = ((18392342566 COS (v) - 11352051608 cos(2t;)+ 
4958070583 cos (3t;) - 1405810666 cos (4t;) + 
234300323 cos(5i;))t;2)/7257600 - ^^Ssoo " 
2cos(4v) +4cos(5'u) - 4cos(6i;) + 2cos(7?;) 

^2,num = -((35142254976 COS (u) - 20245959411 cos (2^;) 
+7950775936 cos (3v) - 1405906674 cos (4?;) 
+234300323 cos (6?;) )w2)/7257600 + ^^^^^^ + 
24cos(4u) - 48 cos (5^;) + 48 cos (6^;) - 24cos(7t;) 

bl,num = ((50246280942 cos (u) - 26679563229 cos (2t;) 
+8977155979 cos (3i7) - 702953337 cos (5w) 
+702905333 cos(6?;))t;2)/3628800 - ^^^fj^^'"' - 
132 cos (4i;) + 264 cos (5u) - 264 cos (6t;) + 132 cos (7?;) 

Knum = -((119523462784 cos (f) - 43206415175 cos (2?;) 
+17954311958 cos (4t;) - 7950775936 cos (5t;) 

+4958070583 cos (6i;))i;2)/7257600 + ^^^^^^+ (33) 
440 cos(4i;) - 880 cos(5?;) + 880 cos(6?;) - 440 cos(7'y) 

bl,num = ((113384696634 COS (i;) - 43206415175 cos (3i;) 
+53359126458 cos (4t;) - 20245959411 cos (5w) 

+11352051608 cos(6i;))t;2)/7257600 - ^^^1209600^"' " 
990 cos (4^;) + 1980 cos (5?;) - 1980 cos (6?;) + 990 cos (7^;) 

^6,n«m = -((56692348317 COS (2t;) - 59761731392 cos (3t;) 
+50246280942 cos (4^;) - 17571127488 cos (5?;) 

+9196171283 cos (6z;))t;2)/3628800 + ^^^604800^^' + 

1584 cos (4i;) - 3168 cos (5i;) + 3168 cos (6i;) - 1584 cos (7t;) 

&7,n«m = (v^ (-7187836062 COS (i;) + 37562934057 COS (2v) 

-36857631107 cos (3t;) + 29909316888 cos (4t;) 
-10358730975 cos (5u) + 5373508799 cos(6t;)))/1814400- 
1848(cos(4w) - 2cos(5u) + 2cos(6u) - cos(7w)) 
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ifl _ 433489274083 _ 152802083671^;^ , 1000430523577^)^ 
"T,1 — 237758976000 2853107712000 291016986624000 

69882256253489?;6 , 257597135900761^)^ 91527043218239^)^° 

1548210368839680000 1532728265151283200000 3384264009454033305600 

ffi _ 28417333297 , 152802083671^)^ _ 1000430523577^)^ , 
"T,2 4953312000 237758976000 24251415552000 

698822562534891)^ 257597135900761^■^ , 91527043218239^;"^ , 

129017530736640000 127727355429273600000 282022000787836108800 

ifi _ 930518896733 _ 1680822920381^^ , 11004735759347^)^ 
"T,3 39626496000 475517952000 48502831104000 

768704818788379i>'^ , 257597135900761i>'^ 91527043218239f i° _ 

258035061473280000 23223155532595200000 51276727415970201600 

ifl _ 176930551859 , 1680822920381d^ _ 11004735759347^)"^ , 
"r,4 — 2971987200 142655385600 14550849331200 

7687048187883791)^ _ 257597135900761i)^ , 91527043218239i)^° , 
77410518441984000 6966946659778560000 15383018224791060480 

;,0 _ 7854755921 _ 1680822920381d^ , 11004735759347d^ 
"T,5 65228800 63402393600 6467044147200 

768704818788379t)S , 257597135900761d^ _ 915270432182391)^° _ 
34404674863104000 3096420737679360000 6836896988796026880 

ifi _ 146031020287 , 16808229203811)^ _ 110047357593471)"^ , 
"T,6 825552000 ^ 39626496000 4041902592000 

7687048187883791)'^ _ 257597135900761d^ , 91527043218239d^° , 
21502921789440000 1935262961049600000 4273060617997516800 "I" • • • 

;,0 _ 577045151693 _ 1680822920381d^ , 110047357593471)-^ 
"T,7 2830464000 33965568000 3464487936000 

7687048187883791)*^ , 257597135900761d^ _ 91527043218239d^" _ 
18431075819520000 1658796823756800000 3662623386855014400 



(34) 
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Method PF-Dl: 
b\,num = (4?; sin^(l) (29030400(2 COS (v) + 2cos(2?;) + 2cos(3t;) 

+2 cos(4f ) + 2 cos(6t;) + 1) sin3(|) + i;(3628800(9 cos(^) 
-19cos(f ) + 2(11 cos(ii^) - 7cos(if^) + cos(i|^))) 
-llv2(65542714cos(|) - 133977068 cos (|^) + 127463860 cos (|^) 
-62185337 cos(^) +21299831 cos(^)))))/14175 
bl,num = (8wsin9(|)(i;(lli;2 (57295722 cos (f) - 50530458 cos (^) 
+72737235 cos (^) + 6776053 cos (^) + 1285617cos(^) 
+21299831 cos(4^)) - 1814400(10 cos (|^) + 68 cos(J^) 
-156cos(|^) + 187cos(ii^) - 119cos(4^) + 9cos(i|^) 
+ cos(i|^))) - 29030400(12 COS (t;) + 12cos(27;) + 12cos(3t;) 
+llcos(4t>) + 2cos(5i;) + 10cos(6t;) + cos(7f ) + 6) sm3(|)))/14175 
b\,num = (4w sin^d) (58060800(66 cos (u) +66 cos(2u) +66 cos(3?;) 
+56cos(4?;) + 20cos(5v) + 46cos(6i;) + 10cos(7?;) + 33) sin^(f ) 
+z;(7257600(50cos(f ) +97cos(^) -267cos(f ) 
+341cos(ii2£) _ 217cos(i|^) - 9cos(if^) + 5cos(i|2^)) 
-11^2(418185576 cos(|) - 101897295 cos(|^) + 429149785 cos(f) 
+213355284 cos (^) + 25712340 cos(|^) 
+21299831(9 cos(ii2l) + cos(if^))))))/14175 

bl,num = (8i; sin9 (|)(t;(lli;2 (469639178 cos (I) + 311586932 cos (1^) 
+333470325 cos (1^) +480049389 cos (^) + 77311321 cos(|^) 
+260311901 cos (ii^) + 63470954 cos (4^)) - 9072000(90 cos (f) 
+36cos(^) - 188cos(|^) + 275cos(ii^^) - 175cos(i|^) 
-47cos(%) + 9cos(i^))) - 145152000(44 cos (t;) +44 cos(2t>) +44 cos(3i;) 
+35cos(4i7) + 18cos(5t;) +26cos(6?;) + 9cos(7'y) + 22) sin3(|)))/14175 
blnum = (4t;sin9(f)(435456000(66cos(t;) + 66cos(2t;) +66cos(3i;) 
+50cos(4t>) + 32cos(5v) + 34cos(6t;) + 16cos(7t;) + 33) sin^d) 
+?;(54432000(80cos(f ) - 23cos(^) - 51cos(f ) 
+110cos(ii^) - 70cos(if^) - 54cos(i|^) + 8cos(^)) 
-llz;2(2105070006cos(f ) + 1324106064 cos (|^) + 1778508400 cos(^) 
+1717441153 cos(^) + 754192017cos(f ) + 920962652 cos (ii^) 
+380999708 cos ( if^ )))))/ 141 75 

binum = (16t; sin9(f)(u(llu2 (809642310 cos(|) + 579296403 cos(^) 

+714780380 cos (^) +616166543 cos (^) +386753499 cos (^) 
+310811926 cos(4^) + 175060939 cos(%)) - 5443200(350 cos(|^) 
-212cos(^) + 12cos(^) + 209cos(^) - 133cos(i|^) 
-261 cos(4^) + 35 cos(^))) - 87091200(132 cos(w) 
+132cos(2t;) + 132cos(3t;) +97cos(47;) + 70cos(5i;) 
+62cos(6f) + 35cos(7z;) + 66)sin3(|)))/14175 

binum = -(8usm9(|)(i;(lli;2(l943141544cos(|) + 1212190059 cos (f) 
+1835374595 cos (f) + 1290288356 cos (^) + 1000981284 cos (f) 
+673851637 cos(ii^) + 426700525 cos(if^)) - 152409600(30 cos (^) 
-21cos(^) + 7cos(f ) + llcos(^) -7cos(if ) 
-23cos(%) + 3cos(i^))) + 304819200(3 sin(f) - 4sin(^) 
+2sm{^) - 4sin(if^) + 3sin(if ))))/14175 
^Lnom = (-4194304^;4cos(|)sm^i(|)) 

(35) 
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1.1 _ 433489274083 _ 152802083671^)^ , 680989543811^;-* 
"T,l — 237758976000 1426553856000 116406794649600 

125177474703917^;^ , 517885739552761i)^ 2572884198423151^;*° _ 

2322315553259520000 306545653030256640000 211516500590877081600000 

Ll _ 28417333297 , 152802083671^;^ _ 1000430523577^;^ , 
"T,2 4953312000 118879488000 8083805184000 

161750007895703f" _ 2419392089643157t)'^ , 69067938626578009^)*° _ 
21502921789440000 6386367771463680000 5875458349746585600000 

tl _ 930518896733 _ 1680822920381^^ , 851496508169t)^ _ 
"T,3 39626496000 237758976000 923863449600 

31098226832101431)" , 17171854137770701u'^ _ 13736401199362907271)*° , 
43005843578880000 ' 4644631106519040000 11750916699493171200000 

tl _ 176930551859 , 1680822920381i;^ _ 7685041522471i;'* , 
"T,4 — 2971987200 71327692800 2078692761600 

373024123233931571)^ _ 1150037153857349d^ , 5553336881578048313i;*° _ 
116115777662976000 69669466597785600 10575825029543854080000 

ul _ 7854755921 _ 1680822920381i;^ , 4465879941727i^ _ 
"T,5 65228800 31701196800 479040307200 

14651758435060069i;g , 5432847035340293i;^ _ 2192163846661534231i;*° , 
17202337431552000 123856829507174400 1566788893265756160000 

ul _ 146031020287 , 1680822920381i;^ _ 8558110979591)-* , 
"T,6 825552000 19813248000 53892034560 

159823310314174791;*' _ 436210741712267i;^ , 798931592780948369i;*° _ 
10751460894720000 5691949885440000 326414352763699200000 

lA _ 577045151693 _ 16808229203811)^ , 1309693001162571)-* 
"T,7 2830464000 16982784000 6928975872000 

492775656906098471;'^ , 335110207212583i;^ _ 7395015266709846197i;*° , 
27646613729280000 3645707304960000 2518053578462822400000 
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Method PF-D2: 

blnum = ^v^sini^(|)(8i;2(ll(-4002729cos(T;) + 2078430 cos(2i;) 

-724279 cos (3v) + 2346178)f2 + 725760(-62 cos(t;) + 59cos(2i;) 
-40cos(3w) + 26cos(4v) -8cos(5i;) + cos(6w) + 35) sin^d)) cos^d) 
+483840sm3(|)(t;(30cos(t;) +30cos(2t;) +30cos(3u) + 13cos(4t>)+ 
18 cos{5v) + 12 cos(6?;) - 5 cos(7w)) + 3(5?; + sm{4v) + sm{7v)))) 
blnum = sini5(|)(-87;2(-9314063t;2 + 572(22949^^- 

120960) cos{v) + 44(1391040 - 428431?;^) cos(2t;) + (9699668i;2 
-46287360) cos (37;) + (27699840 - 7967069?;2) cos(47;) 
-10644480 cos (5t;) + 1108800 cos (6w) + 241920 cos (7t;) - 60480 cos(8?;) 
+35925120) cos^d) - 322560sm3(|)(2D(90cos(t;) + 90cos(2t>) 
+81cos(3?;) + 50cos(4?;) + 49cos(5?;) + 30cos(6t;) - 4cos(7v) - 2cos(8t;)) 
+3(30i; + 8111(31;) + 4sin(4i;) + sin(5i;) + sin(6i;) + 48111(71;) + sin(8i;)))) 
blnum = sini5(|)(10080(8t;(1980cos(i;) + 1961cos(2i;) 

+i680cos(3i;) + 1226cos(4t>) + 1015cos(5t;) + 591cos(6t;) + 50cos(7i;) 
-52cos(8i;) - 3 cos (9^;) + 990) sin^d) - 39cos(^) + 3(13cos( J) 
+32cos( J) - 46cos(|i) +46co8(i|^) - 32cos(%) - 13co8(if ) 
+13cos(4^) + cos(^))) - cosd)((ll(6333473cos(2i;) 
+4157054 cos (3t;) + 3619356 cos (4i;) + 2960054 cos (5i;) 
+724279 cos(6i;))f 2 +4(17018353i;2 +438480) cos(i;) 
+20160(-652cos(2i;) + 445cos(3i;) - 118cos(4i;) - 675cos(5i;) 
+939cos(6t;) + 129cos(7t;) - 104cos(8i;) + 14cos(9i;) + co8(10i;)))i;^ 
+12(2409451f'^ - 110880i;2 + 2520))) 

blnum = sini5d)(3(432759107i;^ - 72576001;^ + 483840) cosd) 

+362880(17cos(f ) -31cos(f ) - 16 cos(f ) + 37cos(f ) 
-37cos(i|^) + 16cos(if^) + 31cos(i^) - 17cos(i|^) - 4cos(^)) 
+i;(i;(ll(103979586cos(f) + 82091598 cos(^) + 62181040 cos(f) 
+36275904 cos (ii^) + 724279(15 co8(if^) + cos(^)))i;2 
+ (1298904167t;2 - 104025600) cos (|^) + 60480(-18 cos(|^) 
-366cos(^) - 1325 cos(|^) +957cos(ii^) + 2370 cos(i|^) 
+718cos(%) - 279cos(4^) + 15 cos(if^) + 8 cos(^))) 
-1935360(1650 cos (i;) + 1612cos(2i;) + 1365co8(3i;) 
+1066cos(4t;) + 819cos(5t;) + 468cos(6i;) + 100cos(7i;) 
-34 cos(8i;) - 6 cos(9i;) + 825) sin^d))) 
H,num = ^i^^sin^^d) (-24(404699571;^ - 10584001;^ + 
70560) cosd) + 60480(-44cos(f ) + 117co8(f ) + 37cos(f ) 
-109cos(f ) + 109cos(i|^) - 37cos(if^) - 117cos(i^) 
+44cos(i|^) + 28cos(^)) + i;(161280(14850cos(i;) + 14318 cos(2i;)+ 
12210 cos (3i;) + 9699 cos (4i;) + 7266 cos (5i;) + 4152 cos (6i;) + 
1125co8(7i;) - 176cos(8i;) - 84co8(9i;) + 7425)sin=^d) 
+i;(-ll(76636238 cos (^) + 62120365 cos (^) + 46335777 cos (^) 
+26765870 cos(ii2^) + 9510034 cos(if^) + 1196642 cos(if^))i;2 
+14(4348800 - 69382489i;2) cos(f ) + 20160(288 cos (f) 
+1506cos(^) + 1900cos(^) - 2112cos(^) - 4320cos(i|^) 
-1913cos(%^) + 339co8(4^) + 60cos(i|^) - 28cos(^))))) 

(37) 
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bl,num = ^t'2sini5(|)(3(521152357t;4 - 16934400t;2 + 
1128960) cos(l) + 241920(7cos(|^) - 33cos(|^) - 20cos(^) 
+41cos(f ) -41 cos(i|^) + 20cos(i|^) + 33cos(i^) - 7(cos(i|^) 
+2 cos(^))) +t>(z;(ll(123210823cos(^) +99307625 cos(^) 
+73495191 cos(l^) + 43125241 cos(ii2^) + 15932285 cos(^) 
+2736371 cos(ifi))?;2 + (1526016833^2 - 73382400) cos (f) 
-40320(774 cos(f) + 1002cos(^) + 1333cos(f ) - 1617cos(^) 
-3258cos(if^) - 1442cos(if^) + 51cos(i|^) + 105cos(i|^) 
-28cos(^))) - 1290240(2970 cos (i;) + 2837cos(2v) + 2445cos(3t;) 
+1938cos(4i;) + 1446cos(5t;) + 831cos(6t;) + 240cos(7t;) - 14cos(8t;) 
-21cos(9f) + 1485)sm3(|))) 

62^^^ = -^v"^ smi5(|)(54(50601001f^ - 1764000?;^ 
+117600) cos(f ) + 1270080(cos(f ) - 9cos(f ) - 8cos(^) 
+14cos(f ) - 14cos(if^) + 8cos(i|^) + 9cos(i|^) - cos(i|^) 
-5cos(^)) + i;(?;(ll(3(71792647cos(f ) + 57772764 cos(l^) 
+42931200 cos(l^) + 25158645 cos(ii2^) + 9410087 cos(if£)) 
+5313226 cos(if^))t>2 + (2670318541i;2 - 112190400) cos(^) 
+423360(- 171 cos ( f ) - 141 cos ( f ) - 221 cos ( f ) + 264 cos ( ) 
+540cos(ifi) +229cos(i|^) + 6cos(i|^) - 21cos(if^)+ 
5cos(^))) - 3386880(1980 cos (w) + 1885cos(2i;) + 1632cos(3i;) 
+1290cos(4w) + 963cos(5i;) + 555cos(6t;) + 162cos(7v) - 4cos(8i;) 
-15cos(9u) + 990) sin^d))) 
^lenam = (2147483648^;^ cos3(|)sin27(f ) 

(38) 
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1.2 _ 433489274083 _ 152802083671^;^ , 1404086671901^;'* 
"T,l — 237758976000 951035904000 194011324416000 

108627551857199^;^ , 3113473234169t)^ 21678565330566029^;^° _ 

1161157776629760000 1621934672117760000 282022000787836108800000 

l2 _ 28417333297 , 152802083671^;^ _ 1000430523577^;^ , 
"r,2 4953312000 79252992000 4041902592000 

3812117933243383t)^ _ 131666706221101^)^ _|_ 766613393985947587t)"' _ 
193526296104960000 133049328572160000 23501833398986342400000 

l2 _ 930518896733 _ 1680822920381;;^ , 67397661839051;;'^ 
"r,3 39626496000 158505984000 32335220736000 

475080967011229691;" , 31127602487128507t)^ _ 59333732949165745199i;i° , 
193526296104960000 ' 1548210368839680000 47003666797972684800000 

l2 _ 176930551859 , 1680822920381;;^ _ 85581 1097959;;-* , 
"r,4 — 2971987200 47551795200 97005662208 

149201016148079837;;^ _ 407769624909121;;^ , 6653867251060213627;;^° _ 
116115777662976000 3225438268416000 742163159967989760000 

l2 _ 7854755921 _ 1680822920381;;^ , 19713857381587;;^ 
"T,5 65228800 21134131200 862272552960 

10823009510563069;;^ , 83749133157903719;;^ _ 189076914789983483663;;^° , 
2867056238592000 206428049178624000 6267155573063024640000 

;,2 _ 146031020287 , 1680822920381;;^ _ 26590548293789;>-* , 
"T,6 825552000 13208832000 673650432000 

224945948304809533;;*^ _ 12256588145611;;^ , 232561853289543390209;;^° _ 
32254382684160000 15672683520000 3916972233164390400000 

1 2 _ 577045151693 _ 1680822920381;)^ _, 1089598285975631)"^ 
"T,7 2830464000 11321856000 2309658624000 

117725260678970569;;^ , 35655584375317913;;^ _ 248038978837339401007;;^° , 
13823306864640000 36862151639040000 3357404771283763200000 

r39 
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Method PF-D3: 
blnum = ^^^;^cos(|)sini8(|)(140734cos(|)i;5 + 12(5357?;4 

-ieSOf 2 + 480) cos{^)v + (4(21791^'' + 6480v'^ - 3600) cos(^) 
+(97229?;'' - 27360^2 + 8640) cos(^) + (32989t;'' - 7200?;^ 
+10800) cos{^) + 80(-8?;(492 cos{v) + 492 cos(2u) + 301 cos{3v) 
+288cos(4i;) + 215cos(5i;) - 86cos(6?;) - 60cos(7t;) + 26cos(8t;) 
+246) sin^d) + 33(16t;2 - 9) cos(i^) + 54(3 - 2v^) cos{^) 
+30(3 - 4i;2) cos(ifi) + 9(8t;2 - 13) cos{^) + 3(9 - 4v^) cos{^)))v 
-61440cos2(f )(2cos(t;) + 2cos(2t>) + 2cos(3z;) + 2cos(4i;) 
+2cos(6t;) + l)sin5(|)) 

blnum = -^^^v^ cos3(§)smi8(|)(84282cos(f)i;5 + 4(20504t;4 
-4200t;2 + 1125) cos{^)v + (12(4774t>'' + 2160^;^ - 975) cos{^) 
+15(4015?;'' - 2112?;2 + 552) cos(^) + (32989?;'' + 4860?;2 
+6480) cos{^) + 60(99(5?;2 - 3) cos{^) + 3(63 - 71?;^) cos(i|^) 
+(60 - llv'^) cos(if^) + 2(-4z;(492 cos(?;) + 492 cos(2t;) 
+286cos(3?;) + 310cos(4i;) + 202cos(5?;) - 95cos(6?;) - 38cos(7?;) 
+19cos(8?;) + 246) sm'^(|) + 3(5z;2 - 17) cos{^) 
-3(f2 - 4) cos(if^))))t> - 46080cos2(|)(2cos(i;) + 2cos(2i;) 
+2cos(3i;) + 2cos(4t;) + 2cos(6?;) + l)sm5(|)) 
blnum = ^^^^^2cos(|)sini8(|)(-61440cos2(|)(22cos(?;) 
+22cos(2?;) + 21cos(3i;) + 16cos(4i;) + 13cos(5i;) + 9cos(6?;) 
+6 cos(7?;) + cos(8?;) + 11) sin^d) + 2i;(604527t;4 - 31000?;2 
+6600) cos(§) + ?;(1090199t;^ - 5920t;2 - 17520) cos{'^) 
+?;((951159?;'^ - 47360?;2 - 7560) cos(|^) + 5(148753?;^ 
-10576?;2 + 3480) cos(f ) + (444741?;" + 48080?;^ - 7200) cos(f ) 
+120(-132cos(ii^^) + 60cos(i|^) + 218cos(i|^) - 36cos(i|^) 
-118cos(if^) + 15cos(^) + 7cos(^)) + ?;(11?;(15863?;2 
+8960) cos{^) + ?;(32989?;2 + 37280) cos(if^)- 
40(8(10601 cos(?;) + 9368008(2?;) + 7755cos(3?;) + 5858cos(4?;) 
+3103cos(5?;) + 538cos(6?;) - 457cos(7?;) - 154cos(8?;) 
+70 cos(9?;) + 14 cos(10?;) + 5412) sin3(|) + ?;(383 cos(if^) 
+103cos(i|^) +4(-26cos(if^) + 2cos(%) +cos(^))))))) 
blnum = -^i^?;3 cos3(|)smi8(f)(-15360cos2(|)(110cos(?;) 
+il0cos(2?;) + 108cos(3?;) + 78cos(4?;) +66cos(5?;) +44cos(6?;) 
+32 cos(7?;) + 2 cos(8?;) + 55) sin^d) + 6?;(252263?;^ - 15900?;^ 
+3720) cos(f ) + ?;(1374791?;'' + 33000?;^ - 30780) cos(f ) 
+?;(5(239129?;'' - 18288?;^ - 2052) cos(f^) + 15(62139?;" 
-4720?;2 + 2040) cos(^) + (553817v" + 79080?;^ - 12600) cos(^) 
+180(-143cos(ii^) + 59cos(i|^) + 237cos(i|^) - 75cos(i|^) 
-104cos(i|^) + 28cos(^) + 2cos(^)) + ?;(165?;(1243?;2 
+760) cos(%) + ?;(32989?;2 + 42360) cos(if^) - 20(8(26818 cos (?;) 
+23302 cos (2?;) + 19472 cos(3?') + 14923 cos(4?;) + 7646cos(5?;) 
+934 cos (6?;) - 1360 cos (7?;) - 200 cos (Su) + 194 cos (9?-) 
+llcos(10?;) + 13530) sin^d) + 3?;(426cos(^) + 5cos(i^) 
_77cos(i|^) + 15cos(^) + cos(f^)))))) 

(40) 
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blnum = ^^i'^cos(|)sini8(^)(-184320cos2(|)(330cos(t;) 

+326cos(2i;) + 298cos(3w) + 250cos(4t;) + 188cos(5t;) + 142cos(6i;) 
+80cos(7i;) + 32cos(8z;) +4cos(9t;) + 165) sin^d) + 480t;(109813i;^ 
-4160t;2 + 384) cos(f ) + 12w(4062663u'' - 134080^2 - 19200) cos(|^) 
+v(12(3490311i;^ - 138560t;2 - 16800) cos(f ) + (32182799^^- 
992640^2 + 192960) cos(^) + (20304031?;^ + 1640640t;2- 
74160) cos(f ) + 720(-627cos(^) + 532cos(if^) + 788cos(i|^) 
-17cos(i5^) - 385cos(i|^) - 176cos(^) + 52cos(^) 
+12cos(4^)) + i;(t;(32989(17cos(^) +cos(if ))?;2 + 176(54469?;2 
+16860) cos{^) + 16(192181?;2 + 97440) cos(%) + 480(352 cos(if^) 
-232cos(i|2£) + 22cos(if^) + 68cos(^) - llcos(^) - 3cos(^))) 
-1920(77156 cos (?;) + 69308 cos (21-) + 57403 cos (3z;) + 42088 cos (4?;) 
+23663 cos (5^;) + 7258cos(6w) - 332cos(7t;) - 962cos(8w) 
-52 cos(9z;) + 152 cos(10t;) + 22 cos(llt;) + 40106) sin^d)))) 
blnum = -M^^3,„,3^.),i^i8(.)(_92160cos2(|)(66cos(^) 
+66cos(2'y) +61cos(3'y) + 50cos(4t;) + 37cos(5?;) + 29cos(6?;) 
+16cos(7u) + 5cos(8t;) + 33) sin^d) + 6t;(889361t;^ - 36900?;2 
+5880) cos(|) + t;(4932719f - 121800t;2 - 42840) cos(|^) 
+?;(3(1410343i;4 - 62640?;^ - 11880) cos(|^) + 48(68101?;^ 
-3465i;2 + 1110) cos(^) + 4(502843w^ + 53175i;2 - 2790) cos(^) 
+3(11(26231?;^ + 109801;^ - 2520) cos(ii^) + (78023t;^ + 46140i;2 
+19560) cos(if^) + 2((4829v'' - 1870^^ + 13890) cos(if^) 
-10(8i;(31262cos(t;) + 27954 cos(2i;) + 23122 cos(3i;) + 17249 cos(4t;) 
+9324cos(5i;) + 2061cos(6f) - 650cos(7f) - 363cos(8'y) + 74cos(9i;) 
+55cos(10i;) + 16236) sin^d) + 3(67t;2 + 99) cos(^) 
+21(32 - 5v^) cos(if^) + (84 - 31^;^) cos(^) 
+15(T;2-6)cos(f^)))))) 

blnum = ^^v^ COS d ) sinis d)(- 122880 cos^d) (462 COS (?;) 
+452cos(27;) +417cos(3t;) +344cos(4t;) + 271 cos(5f ) + 191 cos(6i;) 
+ 118cos(7?;) +45cos(8?;) + 10cos(9f ) + 231) sin^d) 
+2?;(24403159?;'' - 927600?;^ + 75600) cosd) + 2t;(22599203u^- 
730560i;2 - 111600) cos(|^) + t;(6(6473137f'* - 267360i;2 
-17880) cos(^) + 2(14893043t;^ - 320160f ^ + 65880) cos(^) 
+6(3171377^^ + 240960?;^ - 28320) cos(f ) + 720(-264cos(ii^) 
+302cos(i|^) +664cos(if^) + 50cos(i|^) - 348cos(i|^) 
_127cos(4^) + 5cos(^) + 20cos(^)) +i'(f(66(142321t;2 
+36320) cos(ii^) + 2(1641761^2 + 744960) cos(if^) + (7299711;^ 
+262800) cos(%) + (84227t;2 - 61680) cos ( ) + 480(24 cos (i|^) 
+32 cos(^) + 5 cos(^) -5 cos(^))) -640(214989 cos(?;)+ 
192536 cos (2t;) + 159857 cos(3t;) + 116994 cos (4t;) + 66367 cos(5t;) 
+23162 cos(6f) + 1179cos(7v) - 2150cos(8t;) - 278cos(9i;) 
+270cos(10u) + llOcos(llu) + 111232) sin^d)))) 
blnom = (824633720832t;8 cos^ d) sin^o (§)) 

(41) 
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u3 _ 433489274083 _ 152802083671^^ , 2211398968549^;^ 
"T,l 237758976000 713276928000 291016986624000 

335780690096891)" 144902264134913t)^ 18020995400748499t)"' 



145144722078720000 17516894458871808000 14101100039391805440000 

l3 _ 28417333297 , 152802083671^^ _ 1000430523577i;-^ , 
"T,2 4953312000 59439744000 2425141555200 

666660081166011)'^ _ 11606680689206023t)^ , 363627917613911087i;i° _ 
1860829770240000 6386367771463680000 5875458349746585600000 

u3 _ 930518896733 _ 1680822920381^^ , 180183513998459i;^ 
"T,3 39626496000 118879488000 48502831104000 

67733305508864471;^ , 8117004168919561i;^ _ 9618739589821913801d^° , 
12095393506560000 142911726354432000 2350183339898634240000 "r • • • 

i3 _ 176930551859 , 1680822920381i;^ _ 117366928934503i;-* , 
"T,4 2971987200 35663846400 7275424665600 

94400454891172671;^ _ 154456853448146527i;^ , 156768697509684951877i;" _ 
2902894441574400 348347332988928000 3525275009847951360000 

t3 _ 7854755921 _ 1680822920381d^ , 210537222465471)'' 
"T,5 ~ 65228800 15850598400 497464934400 

866895436403651)*^ , 153981351646932977i;^ _ 981460420389037009991;^° , 
8601168715776 95274484236288000 522262964421918720000 • ■ 

ifi _ 146031020287 , 1680822920381i;^ _ 1485385540033871;" , 
"T,6 825552000 9906624000 2020951296000 

770892579458067231;*^ _ 92261723864590011;*^ , 16273137531259548461d^° _ 
4031797835520000 2764661372928000 39169722331643904000 



"T,7 — 



577045151693 _ 1680822920381i;^ , 60974002854799i;'* 
2830464000 8491392000 692897587200 



203359037562761171;^ , 2799280124854146809i;8 _ 449833739846395057357i;i° , 
863956679040000 663518729502720000 839351192820940800000 

(42) 

Method PF - D4: 

binum = 1006632961;^ cos^d ) sin^^d ) ((11858^6 - 1512?;'^ + WObv'^ - 
60) cos(f ) + 9(1540i;6 + 280v^ - 27lv^ + 20) cos(f ) + 60( - cos(f ) - 
6 cos(f ) + 8 cos(f ) - 8 cos(if^) + 6 cos(if^) + cos(i|^) - 3 cos(if^) 
+ cos(^) ) + v{- 3072 cos^d) (18 cos{v) + 18 cos(2t>) + 18 cos(3f ) + 
3 cos(4u) + 14 cos{5v) + 4 cos{6v) - 7 cos(7i;) + 9) sin^d ) - 16i;2(2244 
cos(t;) + 1419 cos(2t;) + 1480 cos(3t;) + 914 cos(4t;) - 961 cos(5?j) - 
535 cos(6?;) + 428 cos(7u) + 104 cos(8^;) - 77 cos(9?;) + 1122) sin^d) 
+ v{llOUv'^ - 3816^2 + 1401) cos(f ) + ?;((4675u^ + ISOu^ + 2358) 
cos(^) + (803v^ + 5028^2 - 5376) cos(^) + 3(44(18 - 13^^) cos(l^) + 
4(283 - 155i;2) cos(lf^) + 2(212f2 - 525) cos(l|^) + (40^2 - 71) cos{^) 
+ (285 - 92v^) cos(if^) + (20^;^ - 71) cos(^) )))) 
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binum = - 100663296^^ cos^d) sm^^d) (2(77077t;6 - 672v'^ - 107Av^ + 
120) cos(|) + 6{23705v^ - 1920^'' - 22v'^ + 40) cos(f ) - 240(4 cos(f ) - 
cos(^) - 2 cos(f ) + 2 cos(if^) + cos(if^) - 4 cos(^) + cos{^) + 2 
cos(^) - cos(^) )+v{- 6144 cos^d) (108 cos(t;) + 108 cos(2i;) + 77 
cos(3?;) + 70 cos(4t;) + 43 cos{5v) + 21 cos(6t;) - 6 cos(7^;) - 13 cos(8t;) 
+ 54) sin^d) - 32i;2(11668 cos{v) + 10196 cos(2?;) + 7776 cos(3t;) + 
2663 cos(4i;) - 1296 cos(5t;) - 1641 cos(6i;) + 164 cos(7t;) + 661 cos(8i;) 
+ 12 cos(9t;) - 107 cos(10t;) + 6732) sm^{^) + 7;(113707t;'^ - 69727;^ + 
7080) cos(l^) + f (67925f'^ + 10140^^ - 8412) cos(^) + w(3(8283w^ + 
4572^2 - 1120) cos(f ) + 11(365^;^ - Uv"^ + 432) cos(^) + 12((742 

- 505i;2) cos(if^) + {57v^ - 124) cos(if£) + (185t;2 - 689) cos{^) + 
(143 - 41i;2) cos(i|^) + 2(86 - 15v^) cos(^) + (lOu^ - 59) cos(^) 
)))) 

H,num = 201326592^4 cos^d) sm^^{^) ((411884i;6 . 10092^^ - 51 15^^ 
+ 660) cosd) + 3(126445i;6 - 7800t;4 + 1598t;2 - 80) cos(f ) + 60( - 14 
cos(f ) - cos(f ) + 13 cos(l^) - 13 cos(i|^) + cos(if^) + 14 cos(i^) 
+ cos(if^) - 10 cos(^) - cos(^) + 3 cos(^) ) + v{ - 3072 cos^d) 
(594 cos{v) + 546 cos(2i;) + 458 cos(3t;) + 355 cos(4t;) + 244 cos{5v) 
+ 108 cos(6v) - 3 cos(7i;) - 40 cos(8t;) - 18 cos(9?;) + 297) sin^d) + 
27;(14895l7;4 - 3768z;2 _^ ^086) cos(f^) + z;(186340z;4 + 201001;^ _ §597) 
cos(^) + v (33 (2648?;'' + 6367;2 - 187) cos{^) + 132(205?;^ + 29t;2 + 81) 
cos(ii^) + ll(365t;i - 450^^ + 1191) cos{^) - 3(731 cos(i|^) + 3298 
cosl^) + 593 cos(i|^) - 1028 cos(^) - 149 cos(^) + 147 cos(^) ) 

- 2?;(8(64294 cos{v) + 54653 cos(2t;) + 38493 cos(3i;) + 15944 cos(4u) - 
1305 cos(5f ) - 4899 cos(6t;) - 676 cos(7f ) + 1690 cos(8t;) + 633 cos(9i;) 

- 230 cos(10t;) - 117 cos(llt;) + 34074) sin^d) + 3v(171 cos(i|^) - 365 
cos{^) - 73 cos{^) + 10(9 cos(^) + cos(^) - cos(^) )))))) 
binum = - 100663296?;^ cos^d) sin^^d) (4(6745751;^ - 23478^^ - 4140i;2 
+ 720) cosd) + 4(610445?;^ - 25050t;'^ + 2328^;^ - 240) cos(|^) + 240( 

- 12 cos(f ) - 9 cos(f ) + 19 cos(f ) - 19 cos(i|^) + 9 cos(if^) + 11 
cos(i|^) + cos(i|^) - 7 cos(^) - 5 cos(^) + 3 cos(^) + cos(^) 
) + v{- 6144 cos^d) (1947 cos{v) + 1780 cos(2t>) + 1532 cos{3v) + 
1153 cos(4t>) + 812 cos(5u) + 376 cos(6z;) + 35 cos(7f ) - 80 cos(8t;) - 64 
cos(97;) - 11 cos(10?;) + 990) sin^d) + i;(1934185f'' - 35376?;2 - 1440) 
cos(f ) + t;(1264615v'' + 92640^^ - 22932) cos(f ) + i;( (642895?;'' + 
114648^2 - 39588) cos(|^) + 11(21955?;'' + 2712t;2 + 5508) cos{^) 
+ 2( - 16?;(210376 cos(?;) + 176605 cos(2?;) + 122822 cos(3?;) + 57558 
cos(4?;) + 5448 cos(5?;) - 9243 cos(6?;) - 2316 cos(7?;) + 2721 cos(8?;) 
+ 1892 cos(9w) - 50 003(101;) - 362 cos(ll?;) - 61 cos(12?;) + 109790) 
sin^d) + (30745?;4 - 5910?;^ _^ 33352) cos{^) + (4015?;*^ - 1242?;^ - 
8928) cos(if^) + 6(5(87?;2 - 608) cos{^) + (261z;2 - 1220) cos(^) + 
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5(130 - 17v'^) cos(^) + (499 - 75^2) cos(^) + 3{5v^ - 39) cos(^) 
+ (5t;2 - 41) cos(^) ))))) 

bl^^^ = 100663296t;4 cos^d) sin^^d) (5(1193060^;^ - 41904^^ - 5997?;2 
+'l020) cos(|) + 60( - 14 cos(f ) - 126 cos(f ) + 3 cos(^) + 68 cos(|^) 

- 68 cos{^) - 4 cos(i|^) + 113 cos(i|^) + 9 cos{^) - 63 cos(^) - 22 
cos(^) + 5 cos(^) + 13 cos(^) + cos(^) ) + vi - 3072 cos^d) 
(8642 cos(t;) + 8000 cos(2t;) + 6715 cos(3?;) + 5297 cos(4d) + 3528 
cos(5z;) + 1774 cos(6f ) + 335 cos(7t;) - 291 cos(8t;) - 234 cos(9t;) - 84 
cos(10i;) - 5 cos(ll?;) + 4438) sm^{^) + 6t;(900240t;^ - 34950^;^ + 1097) 
cos(f ) + 2u(2142800t;^ - 21390i;2 + 4761) cos(f ) + v((2856205v'^ + 
170220i;2 - 81459) cos(^) + 3(511335^^ + 68340^;^ - 6064) cos(f ) 
+ 44(14405?;^ + 1962i;2 + 2133) cos{^) + 3(37216 cos{^) - 1456 
cos(i|^) - 26299 cos(i^) - 9051 cos{^) + 3561 cos(^) + 2870 
cos(^) + 257 cos(^) - 503 cos(^) - 35 cos(^) ) + t;(t;(20(9515?;2 - 
318) cos(i|^) + 9(4235^;^ - 684) cos(^) + 5(803t;2 + 2052) cos(^) + 
12(404 cos{^) - 84 cos(^) - 100 cos(^) - 12 cos(^) + 15 cos(^) 
+ cos(^) )) - 16(923815 cos{v) + 780435 cos(2t;) + 544332 cos(3?;) + 
265090 cos(4u) + 54870 cos(5v) - 23058 cos{Qv) - 10352 cos(7?;) + 8726 
cos(8t;) + 6750 cos(9t;) + 651 cos(lO'y) - 1000 cos(ll'y) - 444 cos(12t;) 

- 25 cos(13?;) + 488520) sin-'^d) )))) 

bl^^^ = - 100663296t;^ cos^d) sin^^d) (4(2363372z;6 - 84162t;4 - 
11175?;2 + 1860) cosd) + 24(355960^^ - 12795^;"^ + i76v^ - 80) cos{^) 
+ 240( - 39 cos(l^) - 3 cos(^) + 23 cos{^) - 23 cos(i|^) + 2 cos{^) 
+ 35 cos{^) + 6 cos(i|^) - 21 cos(^) - 10 cos(^) + 2 cos(^) + 4 
cos(^) + cos(^) ) + vi - 6144 cos^d) (6889 cos{v) + 6340 cos(2t;) 
+ 5369 cos{3v) + 4201 cos{4v) + 2829 cos(5f) + 1439 cos(6t;) + 331 
cos(7i;) - 177 008(8^;) - 180 cos(9t;) - 69 cos(lOu) - 10 cos(ll?;) + 3530) 
sin^d) + 4t;(1700006i;'' - 14241^^ - 1503) cos(f ) + 4w(1147190i;4 + 
60825^2 - 25641) cos(^) + t>(12(210419t;^ + 24831f2 - 1957) cos{^) 
+ 44(24817t;4 + 3201t;2 + 3159) cos(ii^) + (356345f'^ + 7080t;2 + 
156144) cos(i|^) - 12(262 cos{^) + 8737 cos{^) + 3825 cos{^) - 
2(498 cos(^) + 523 cos(^) + 67( cos(^) - cos(^) )) + 35 cos(^) 
) + t;(t;((82599t;2 - 6408) cos(if^) + 15(803^7^ + 864) cos(i|^^) + (803?;^ 
+ 7584) cos(i|^) - 12(81 cos(^) + 135 cos(^) + 23 cos(^) - 15 
cos(^) - 4 cos(^) )) - 32(731642 cos{v) + 616515 cos(2i;) + 431394 
cos(3v) + 217139 cos(4?;) + 54930 cos(5?;) - 9918 cos(6t;) - 7072 cos(7i;) 
+ 5350 cos(8?;) + 4914 cos(9w) + 747 cos(10i;) - 614 cos(lli;) - 339 
cos(12i;) - 50 cos(13t;) + 386010) sin^d) )))) 

54^^^^ = 201326592^"^ cos^d) sin^^d) (44(14828?;"^ + 1779^2 + 1971) 
cos(ii^) + 48(491i;2 + 42) sind) v + 2(2652 - 18397?;^) sm(f ) v - 
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2(46961^2 + 1572) sm(f ) v - 2(84467i;2 + 996) sin(^) v + 92(132 - 
1315^2) sin(^) V + 12(6461^2 . 1752) sm{^) v + 12(10789f2 - 2244) 
sm(if^) V + 16(2921z;2 + 987) sin(if^) v + 112(219 - 122v'^) sm{^) v 
+ 16(570 - 1019t;2) sm(i|^) v + 8(233^;^ - 1146) sm(^) v + 4(1075t;2 

- 2004) sin(^) v + 10(25t;2 + 12) sm(^) v + 2(372 - 131?;^) sin(^) 
V + 30(12 - 5f2) sm(^) v + (5494852^^ - 200688i;^ - 24198i;2 + 4200) 
cos(|) + (4966940u6 - 164820?;'' + 6873?;2 - 1500) cos(f ) + (3957316?;^ 

- 37644t;^ - 156391;^ - 3900) cos{^) + (2678500t;'' + 135480?;^ - 374131;^ 

- 2340) cos(^) + 2(742786^6 + 87156f'' - 15261t;2 + 2220) cos{^) + 
(217855?;6 + 73201;"^ + 91506i;2 - 4440) cos(i|^) + 4(13079?;^ - 318?;^ 

- 2646^2 + 540) cos(i|^) + 2(4180?;^ + 2820v^ - 253471;^ + 1740) 
cos(i|^) + (737t;6 + 4776t;4 - 28968i;2 + 960) cos(if^) - 6(54?;^ - 863i;2 
+ 380) cos(^) - 24(45f4 - 347^;^ + 80) cos(^) + 9( - 8v^ + 35?;^ + 
60) cos(^) + 3(20^4 - 197i;2 + 140) cos(^) + 9(4?;^ - 35?;^ + 20) 
cos(f^) ) 

bienom = ' (316659348799488?;io cos^o (|) sin^o (|)) 

1.4 _ 433489274083 _ 152802083671t)^ , 77626 18237^)'* 7881601960439^)^ 
"T,l — 237758976000 570621542400 1119296102400 " 14744860655616000 

27284304529514897^)^ _ 1799866965050155021t;i° _ 
613091306060513280000 282022000787836108800000 

t4 _ 28417333297 i 152802083671t)^ _ 1000430523577^)'^ , 604487352966331^)^ 
"T,2 ' 4953312000 47551795200 1616761036800 11058645491712000 

75851624289432059^)^ , 646544241473169703^^° 
25545471085854720000 7833944466328780800000 " • ' ' 
l4 _ 930518896733 _ 1680822920381t)^ , 2349705253321^^ 23296554826706981^;^ 
"T,3 39626196000 95103590400 104190259200 " 22117290983424000 

I 585913207 11987:-{37'»s _ 14407929187812 12()8197i)"' , 

488908537528320000 15667888932657561600000 ' • • • 
1 4 _ _ 176930551859 _, 1680822920381t)^ _ 74576374036553f'^ , 95021198062331d^ 
"T,4: 2971987200 28531077120 2910169866240 14455745740800 

_ 1557322122991096859t)^ 1918393406379510690887»"' _ 

1393389331955712000 14101100039391805440000 
v4 _ 7854755921 _ 1680822920381^)^ , 7297045929049^)"^ _ 20692039318485463^>'^ 
"T,5 65228800 12680478720 107784069120 982990710374400 

I 5526609376838648143;)^ 43203895792 15898805647^ " , 

1238568295071744000 " 6267155573063024640000 "r • • • 
1,4 _ _ 146031020287 , 1680822920381^^ _ 1177252560689^^ , 74732313119187721^^ 
"T,6 825552000 7925299200 9980006400 1843107581952000 

_ 37277993693096489391)'^ , 6574125730067577575911t)"' _ 

387052592209920000 ' 3916972233164390400000 
i4 _ 577045151693 _ 1680822920381^)^ , 122443866047771)-^ 264047572988562471)6 
"TJ 2830464000 6793113600 ' 86612198400 " 526602166272000 

_l_ 81877808195686092431)^ 7493224716658621457999d^° , 

663518729502720000 " 3357404771283763200000 • • • 

Method PF - D5: 

bi,num = - 7247757312^;^ cos^^) sin^^d ) ( - 122880 cos^(|) (2 cos{v) + 
2 cos(2v) + 2 cos(3f) + 2 cos(4t;) + 2 cos(6z;) + 1) sin^(|) + 70i;(396f6 
+ 72v'^ - 95i;2 + 20) cos(|) + 30f(616t;6 - 504v^ + 335v^ - 20) cos{^) 
+ v{60{3v'^{Uv^ + 52v'^ + 35) - 80) cos(f^) + 30(66?;^ + 408t;^ - 
979^2 + 220) cos(f ) + 100(6 cos(f) - 88 cos{^) + 60 cos{^) + 
33 cos{^) - 51 cos(i|^) + 6 cos(if^) + 12 cos(^) - 4 cos(^) ) + 
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v{ - 2560 cos^d) (84 cos{v) + 84 cos(2f;) - 13 cos(3f) + 84 cos(4w) + 
cos(5t;) - 115 cos(6i;) + 31 cos(8i;) + 42) sin^d) - 16z;2(641 cos{v) + 
2520 cos(2i;) - 93 cos(377) - 5564 cos(4?;) + 13 cos(5t;) + 3882 cos(6i;) 

- cos(7i;) - 1292 cos(8t;) + 174 cos(10t;) + 1260) sin^d) + bvHUv"^ 

- 2808^2 + 3966) cos(f ) + 88(44 - 9^;^) cos(ii^) + 24(73^^ - 212) 
cos(^) - 3(24?;2 + 121) cos(i|^) + 3(735 - ISAv'^) cos{^) + 6{20v'^ 

- 71) cos(i|^) + 12{6v^ - 29) cos(^) + 4(29 - Gv"^) cos(^) )))) 
binum = 2415919104t;5 cos^d) sin^^d) ( - 737280 cos4d)(12 cos{v) 
+ 12 cos(2t;) + 12 cos(3i;) + 7 cos(4t;) + 10 cos(5t;) + 2 cos(6?;) + 
5 cos(7v) + 6)sin^d) + 180u(4620?;6 . Q72v^ + 13^^ + 120) cosd) 
+ 45w(14256f6 + 23Uv* + 1987t;2 - 1300) cos(|^) + t;(45(8184D^ 
+ 1608i;^ - 6581?;2 + 1420) cos(f^) + 45(3256^7^ - 588^7"^ + 4583?;^ 

- 420)cos(^) - 300(123 cos(f ) - 154 cos{^) + 130 cos{^) - 96 
cos(i|^) - 84cos(if^) + 171 cos(i|^) - 33 cos(^) - 49 cos(^) + 19 
cos(^) ) + v{- 7680 cos^d) (1008 cos(t;) + 497 cos(2z;) + 814 cos(3u) 
+ 140 cos(4t;) - 191 cos{5v) - 231 cos(6f) - 364 cos(7?;) + 62 cos(8t;) 
+ 137 cos(9t;) + 504) sin^d) - 16i;2(79446 cos{v) + 26700 cos(2t;) - 
54673 cos(3?;) - 35424 cos(4?;) - 6282 cos(5?;) + 23632 cos(6?;) + 23514 
cos(7f) - 7782 cos(8t') - 11000 cos(9t;) + 1044 cos(107;) + 1755 cos(llz;) 
+ 14510) sin^d) + 15?;(3(792t;4 - 652t;2 + 1075) cos(|^) + 22(12^;^ - 
102i;2 - 47) cos(ii^) + 2(582t;2 + 1891) cos{^) + 32(73t;2 - 342) 
cos(i|^) + 24(43 - 36?;2) cos{^) + 3(2131 - SOOv'^) cos(if^) + (364?;2 

- 1887) cos(^) + (132?;2 - 1151) cos(^) + (461 - 60i;2) cos(^) )))) 

= - 4831838208^5 cos^d) sin^-'^d) ( - 368640 cos^d) (66 cos(?;) 
+ 66 cos(2f) + 56 cos(3t;) + 56 cos(4w) + 30 cos(5t;) + 36 003(6^;) 
+ 10 cos(7'y) + 10 cos(8u) + 33) sin^d) + 120?;(17424?;6 - 765^;^ + 
716v^ - 155) cosd) + 45u(37026i;6 + 92v^ - 5159?;2 + 1100) cos{^) + 
i;(15(70422t;6 + 10188i;^ + 9929i;2 - 4220) cos(^) + 300(161 cos(^) 
+ 33 cos(|i) - 319 cos(ii^) + 275 cos(i|^) + 81cos(if^) - 123 cos(i^) 
+ 20 cos(i|^) - 72 cos(^) + 35( cos(^) + cos(^) ) - 18 cos(^) ) 
+ v{ - 7680 cos^d) (2233 cos(u) + 2261 cos(2t;) + 1114 cos(3?;) + 831 
cos(4w) - 199 cos(5t;) - 879 cos(6f) - 369 cos(7t;) - 102 cos(87;) + 137 
cos(9z;) + 121 cos(10z;) + 1386) sin^d) - 16?;2(106239 cos{v) + 46080 
cos(2z;) - 45422 cos(3z;) - 97266 cos(4w) - 7053 cos(5v) + 48138 cos(6u) 
+ 23511 cos(7t;) - 2538 cos(8i;) - 10990 cos(9t;) - 5154 cos(10t;) + 1755 
cos(llt;) + 1270 cos(12z;) + 93890) sin^d) + 15v{{36v'^ (dGSv"^ + 57) 

- 15475) cos(^) + (12672z;4 - 8132^;^ + 20199) cos(|^) + 11(270?;"^ - 
138^2 + 1105) cos(ii^) + (330t;^ + 3138i;2 - 18373) cos(if^) + 3(284?;2 

- 921) cos(if^) + 3(453 - 68^2) cos(^) + 2(1081 - 267^;^) cos(^) 
+ 6(408 - 41t;2) cos(^) + (228?;^ - 1609) cos(^) + 5(12^;^ . 131) 
cos(^) + 4(93 - lOv^) cos(^) )))) 
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binum = 12079595520?;5 cos^d) sm^^{^) ( - 737280 cos^(|) (44 cos{v) 
+'42 cos(2t;) + 40 cos(3w) + 31 cos(4i;) + 28 cos(5w) + 16 cos(6?;) + 
13 cos(7?;) + 4 cos(8t;) + 2 cos(9t>) + 22) sin'^d) + 24t>((4356t;2(25t;2 

- 1) - 1043)t;2 + 530) cos(|) + 3u(716320t;^ + 20196?;'' + 11847^^ - 
11460) cos(f ) + u(9(ll(14520?;^ + 948^2 - 2345)^;^ + 4020) cos(f ) 

- 60(101 cos(^) + 501 cos(f ) - 220 cos(^) + 104cos(^) - 630 
cos(i|H) _ 82 cos{^) + 649 cos{^) - 27 cos(^) - 91 cos(^) - 67 
cos(^) - 30cos(^) + 34 cos(^) )+v{- 1536 cos^d) (16324 cos(i;) 
+ 12089 cos(2t>) + 9812 cos(3z;) + 3800 cos{4v) - 867 cos(5i;) - 3291 
cos(6v) - 3184 cos(7?;) - 484 cos(8i;) + 561 cos(9t;) + 484 cos(10z;) + 
214 cos(llt;) + 8102) sin^d) + 39i;(132(150i;'' + v'^) + 4003) cos(^) + 
v{ - 16t;(160098 cos{v) + 48324 cos(2t;) - 56061 cos(3?;) - 67128 cos{4v) 

- 17958 cos(5v) + 33644 cos(6t;) + 30234 cos(7v) - 534 cos(8?;) - 10044 
cos(9w) - 4320 cos(10t;) + 15 cos(llt;) + 1016 cos(12i;) + 396 cos(13f ) + 
85598) sin^d) + 3(44(2490^2 - m)v^ + 45675) cos{^) + 132(810^^ 

- 276v^ + 509) cos(ii^) + 3(4(1980?;^ + 1716?;^ - 4475) cos(i|^) + 
(880^4 + 8816?;2 - 53214) cos(^) + 2(2861 - 648v^) cos(i^) + (21575 

- 2808i;2) cos{^) + (1959 - 3Uv^) cos(^) + (300?;^ - 2129) cos(^) 
+ (300t;2 - 2501) cos(^) + 2(10^2 - 177) cos(^) + 2(307 - 30v^) 
cos(f^) ))))) 

4,num = - 12079595520^5 cos^d) sin^^d) ( - 368640 cos^d) (2 cos{v) 
+ 1)2(22 cos(i;) + 22 008(2?;) + 20 cos(3v) + 18 cos{4v) + 10 cos{5v) + 
12 cos{6v) + 4 cos(7i;) + 2 cos(8f) + 11) sin'^d) + 6t;(t;2(8580(110i;2 

- 3)?;2 + 1993) - 460)cosd) + 18t;(llt;2 (360(66?;'' + v^) - 827) - 20) 
cos(l^) + ?;(6((13200?;2(41?;2 + 2) - 20807)?;^ + 980) cos{^) + 60(117 
cos(^) - 459 cos(f ) - 440 cos{^) + 640 cos(^) + 783 cos{^) - 
61 cos(i|2^) - 658 cos(i|^) - 336 cos(^) + 217 cos(^) + 169 cos(^) + 
18 cos(^) - 20 cos(^) - 16 cos(^) ) + v{- 1536 cos^d) (34140 cos(t;) 
+ 28661 cos(2?;) + 19811 cos(3?;) + 9491 cos(4?;) - 979 cos(5?;) - 6567 
cos(6?;) - 5641 cos(7?;) - 1793 003(8?;) + 841 cos(9?;) + 1009 cos(10?;) 
+ 428 cos(ll?;) + 95 cos(12?;) + 18514) sin^d) + 9?;(220(933?;2 + 
10)?;2 + 9447) cos(^) + 3?;(484(585?;2 - 64)?;2 + 111477) cos(f) + v{ - 
16?;(304575 cos(?;) + 111318 cos(2?;) - 82171 cos(3?;) - 124536 cos(4?;) - 
30912 cos(5?;) + 53088 cos(6?;) + 50880 cos(7?;) + 6120 cos(8?;) - 15695 
cos(9?;) - 9480 cos(10?;) - 669 cos(ll?;) + 1476 cos(12?;) + 792 cos(13?;) 
+ 162 cos(14?;) + 192432) sin-'^d) + 264(1170?;^ - 204?;^ + 535) cos{^) 
+ 24(3630?;" + 1092?;^ - 7165) cos(if^) + 3((5940?;' + 12544?;^ - 66585) 
cos(^) + (660?;^ + 816?;2 - 5245) cos{^) + 2(15655 - 2208?;^) cos(^) 
+ 4(3141 - 376?;2) cos(^) + (624?;^ - 5363) cos(^) + (456?;2 - 4151) 
cos(^) + 2{32v^ - 303) cos(^) + 16(28 - 3?;^) cos(^) + 4(65 - 6?;^) 
cos(^) ))))) 
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b^nwrn = 2415919104?;5 cos^d) sm^^d) ( - 737280 cos^d) (782 cos{v) 
+ 756 cos(2t>) + 681 cos(3w) + 596 cos(4t;) + 446 cos(5t;) + 346 
cos(6z;) + 196 cos(7w) + 111 cos(8i') + 36 cos(9t;) + 10 cos(lOf) + 
cos(llf) + 395) sin^(l) + 30u(1475232i;6 - 29700^;"^ + 2389t;2 - 2620) 
cos(|) + 15?;(2477376i;6 + 16060i;^ - 104631^2 + 10740) cos(f ) + 
?;(45(578688z;6 + 24420i;^ - 21091;^ - 3460) cos(|^) + 90(168256^^ + 
2904i;4 - 93t;2 + 1500) cos(^) - 300(366 cos(^) + 1240 cos{^) - 1456 
cos(i|^) - 1008 cos(i|^) + 376 cos(i^) + 568 cos{^) + 624 cos(^) 

- 226cos(^) - 306 cos(^) - 15 cos(^) + 17 cos(^) + 23 cos{^) 
+ 2, cos{^) ) + v{ - 7680 cos^d ) (53563 cos(v) + 45974 cos(2w) + 
30514 cos(3z;) + 15658 cos(4z;) - 559 cos(5i;) - 9759 cos{Qv) - 7958 
cos{lv) - 3102 cos{Sv) + 869 cos(9t;) + 1533 cos(10?;) + 666 cos(ll?;) 
+ 190 cos(127;) + 17 cos(13?;) + 29210) sin^d) - 16t;2(2248539 cos{v) 
+ 884628 cos(2z;) - 465296 cos(3z;) - 896802 008(4-?;) - 222138 cos(5f) 
+ 364588 cos(6i;) + 334950 cos(7t;) + 63936 cos(8t;) - 93282 cos{Qv) 

- 73236 cos(lOu) - 9234 cos(llw) + 9682 cos(12z;) + 5604 cos(13t;) + 
1620 cos(14v) + 137 cos(15i;) + 1505344) sin^d) + 30t;(88(2754'i;2 . 
2m)v'^ + 86133) cos(f ) + 15i;(16(11814?;^ - 1449?;^ + 5183)cos(^) 
+ 8(7227t;^ + 1638f2 - 13363) cos(if^) + 8(1683t;^ + 1715t;2 - 9429) 
cos(if^) + 8(29717-^ + 315i;2 _ 1181) cos{^) + 8(33f'^ - 585f2 + 4033) 
cos(i|^) + 8(2988 - 365^2) cos(^) + 10(72^2 - 611) cos(^) + 18(40^;^ 

- 379) cos(^) + (76^;2 - 681) cos(^) + (343 - 36^;^) cos(^) + (385 

- 36i;2) cos(^) + (45 - 4z;2) cos(^) )))) 

= - 4831838208i;5 cos^d) sin^^d) ( - 737280 cos^d) (457 cos(t;) 
+ 436 cos{2v) + 406 cos(3t;) + 331 cos(4t;) + 281 cos{bv) + 181 cos{<6v) 
+ 131 cos(7?;) + 56 cos(8?;) + 26 008(9?;) + 5 003(101;) + cos(lli;) + 
230) sin^d) + 30t;((3960t;2(216?;2 - 5) - 8761)t;2 + 1780) cosd) + 
15u(1435984?;6 + 22000i/ - 3297f2 - 10980) cos(|^) + w(45(337392t;6 
+ 11440t;'i - 27653i;2 + 4460) cos{^) - 300(219 cos(^) + 603 cos(|^) - 
300 cos(^) + 28 cos(i|^) - 756cos(i|^) - 148 cos{^) + 776 cos{^) 
+ 96 cos(^) - 83 cos(^) - 115 cos(^) - 63 cos(^) + 31 cos(^) + 7 
cos(^) + 3cos(^)) + v{ - 7680 cos^d) (32205 cos(w) + 25293 cos(2u) 
+ 19114 cos(3i;) + 8361 cos(4i;) - 241 cos(5z;) - 4841 cos(6i;) - 5127 
cos(7v) - 1558 cos(8t;) + 452 cos(9?;) + 766 cos(lOw) + 452 cos(llz;) 
+ 95 cos(12i;) + 17 cos(13?;) + 16488) sin^d) + 45z;(1408(141z;'' + 
t;2) + 24571) cos(^) + 15i;(352(820t;2 - 49)1;^ + 65325) cos(|^) + v{ - 
16i;(1391199 cos{v) + 496638 cos(27;) - 298296 cos(3t;) - 433962 cos(4w) 

- 149583 cos(5?;) + 183078 cos(6?;) + 208215 cos(7i;) + 31836 cos(8i;) 

- 54047 cos(9t;) - 36606 cos(lOf) - 7569 cos(IId) + 4842 cos(12t;) + 
3624 cos(13i;) + 810 cos(14t;) + 137 cos(15f) + 794124) sin^d) + 
90(189751;-^ - 2632^;^ + 4838) cos (ii^^) + 30(18117?;'' + 2184?;^ - 10934) 
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cos(i|^) + 15( - 2372 cos{^) + 25024 cos(i|^) + 6864cos(^) - 925 
cos(^) - 3401 cos(^) - 1173 cos(^) + 509 cos(^) + 125 cos(^) 

+ 45 cos(^) ) + 30((36(22i>2 + 7) cos(i|^) + 9{llv^ - 168) cos(i|^) 
+ (Uv^ - 488) cos(^) + 48 cos(^) + 192 cos(^) + 56cos(^) - 
2(12 cos(^) + 3 cos(^) + cos(^) ))v'^ + (4488?;'^ + 5348^;^ - 33558) 
cos(i|^) ))))) 

blnom = (-151996487423754240^;i2 cos^^ (|) sin^^ (§)) 

l5 _ 433489274083 _ 152802083671^;^ , 1017850218043^;'* 355108221471443^;^ 
"T,l — 237758976000 475517952000 194011324416000 " 331759364751360000 

131687699860605701^* 9701300523880595811)1" 



1021818843434188800000 47003666797972684800000 

„2 innn/i'^nc;9cic;77„4 2072463900685193«;'5 



i5 _ 28417333297 , 152802083671-»^ 1000430523577<;-^ , 



T,2 4953312000 ^ 39626496000 1154829312000 ^ 27646613729280000 

41477308145052191;* _|_ 25097056509899527ii"^ 



886995523814400000 ' 559567461880627200000 
i5 _ 930518896733 _ 1680822920381^;^ , 270959894639173i;^ 97479391651340473i;g 
"T,3 39626496000 79252992000 32335220736000 " 55293227458560000 

I 1103582448711358933t)* _ 1354275045640832303511?^° , 
5160701229465600000 7833944466328780800000 
^5 _ 176930551859 , 1680822920381-i-^ 180938567211709t)'* , 192417404089068163t)'^ 
"T,4 " 2971987200 23775897600 " 4850283110400 16587968237568000 

13040661300795157i>* _|_ 7536908007008312598671)^° 



5582489310720000 ' 2350183339898634240000 
i5 _ 7854755921 1680822920381;;^ , 60974002854799i;^ 31108033258478857i;^ 
"T,5 65228800 ' 10567065600 615908966400 " 819158925312000 

I 206147997444225374991;* _ 283489566000723918761t)i° , 

2064280491786240000 149217989834833920000 ~r • • • 

146031020287 , 1680822920381d^ 232891275659849i;-* , 3404020481527719231;^ 



;,5 _ _ 146031020287 , 1680822920381d^ _ 232891275659849i;-* _, 

^T,6 825552000 6604416000 1347300864000 4607768954880000 

_ 1791871329414738589D* , 3257163476890690029371d"' _ 

80635956710400000 652828705527398400000 
^5 _ 577045151693 1680822920381?;^ , 4787707284317331)-* _ 3618614330422788731;*^ 
"T,7 2830464000 " 5660928000 2309658624000 3949516247040000 

I 317654346452495203991;* 3797117763219719452879d*° , 

1105864549171200000 ' 559567461880627200000 "r • • • 

Method PF - D6: 

bi,num = - 3478923509760i;*5 cos^^d) sin^^^D {v{ - 12(1459?;'' - 970^2 

+'60) cos(t;) + 4(7379?;'' + 2805v'^ - 1890) cos(2t;) + 180(31 cos(3f) + 
44 cos(4t;) - 78 008(5?;) - 10 cos(6i;) + 88 cos(7?;) - 9 cos(8?;) - 46 cos(9i;) 
+ 4 cos(10?;) + 9 cos(ll?;) ) + ?;(?;( (461 12?;^ - 46950) cos(3?;) + 30(764 
cos(4?;) + 2470 cos(5f ) - 1370 cos(6?;) - 1908 cos(7?;) + 707 cos(8?;) + 
726 cos(9?;) - 124 cos(10?;) - 111 cos(ll?;) ) + ?;(2t;(31457280t;(7 cos(?;) 

- 1) sin^d) cos^^d) - 33436 cos(4?;) - 24673 cos(5?;) + 23310 cos(6?;) 
+ 15023 cos(7?;) - 7754 cos(8?;) - 4901 cos(9t;) + 1044 008(10?;) + 669 
cos(ll?;) ) + 19950 sin(?>) - 5752 sm(2?;) - 58277 sm(3?') + 59598 sm(4?;) 
+ 77003 sin(5?;) - 58792 sin(6?;) - 51447 sin(7?') + 23233 sin(8?;) + 17843 
8111(9?;) - 3480 sin(10?') - 2552 8in(ll?;) )) - 400(24 cos(?;) + 24 cos(2?;) 

- 101 008(3?;) + 108 co8(4?;) - 67 008(5?;) - 173 co8(6?;) + 48 cos(7?;) + 
59 co8(8?;) + 12) sm^{v) )) + 2(6278?'5 - 5265?;''^ + 1170?- + 840 8m(2?;) 

- 210 sin(3?;) - 1260 sin(4?;) + 840 8in(5?;) + 840 8in(6?;) - 1260 8in(7?;) 

- 210 sin(8?;) + 840 sin(9?;) - 210 sin(ll?;) )) 
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^Xnurn = 3478923509760v6 cosiO(|) sin^^^f) (f (8(41029^;'' - 6060^^ - 
2340) cos(w) + (216658?;^ - 231600^2 + 29520) cos(2v) + 720(8 cos{Zv) 

- 73 cos{iv) + 46 cos(5i;) + 20 cos(6u) - 26 cos(7u) + 58 cos(8t;) - 8 
cos(9?;) - 53 cos(10?;) + 6 cos(llt;) + 13 cos(12i;) ) + ?;(z;(4(47520 - 
64489^;2) cos(3u) + 240(889 cos(4v) - 472 cos(5u) + 142 cos(6t;) - 427 
cos(7z;) - 682 cos(8f) + 393 cos(9t;) + 383 cos(10t;) - 84 cos(llt;) - 71 
cos(12t;) ) + i;(2t;(15728640t;(42 cos(2t>) + 47) sin^d) cos^^d) - 53557 
cos(4v) - 23336 cos(5?;) - 10180 cos{6v) + 62776 cos(7t;) + 31222 cos(8t;) 

- 28732 cos(9u) - 15557 cos(lOv) + 4554 cos(llu) + 2637 cos(12u) ) - 
3(4524 sm{v) + 89649 sm{2v) - 109162 sm(3t;) - 66107 sm{4v) + 12508 
sin{5v) - 17082 sm{6v) + 55388 sm{7v) + 44508 sm{8v) - 32396 sin(9t;) 

- 22567 sin(lOf) + 5770 sin(llt;) + 3929 sm(12?;) ))) - 2400(48 cos(t;) 

- 83 cos(2w) + 134 cos(3z;) - 160 cos(4w) - 76 cos(5t;) - 19 cos(6?;) - 
100 cos(7t;) + 46 cos(8f) + 54 cos(9t;) + 24) sin^{v) )) - 4(28927?/'' - 
18240t;3 + 1080?; - 2520 sin(?;) + 630 sin(2?;) + 1260 sin(3?;) - 1890 
sin(4i;) + 1260 8111(5?;) + 1260 sm(6i;) - 1890 sin(7?;) + 1260 sin(8?;) + 
630 sm(9?;) - 1890 sm(10?;) + 630 sin(12?;) )) 

blnum = - 3478923509760i;6 cos^°(|) sin^^^i) (4(221284?;^ - 37695?;^ 

- 9090?; - 1575 sin(?;) + 5670 sin(2?;) + 2520 sin(3?;) - 10080 sin(4?;) 
+ 4095 sm(5?;) + 4095 sin(6?;) - 10080 sin(7?;) + 2520 sin(8?;) + 4095 
sin(9?;) - 1575 sm(10?;) + 2520 sin(ll?;) - 1575 sm(13?;) ) + w( - 2(29854?;^ 
+ 91035?;2 - 24390) cos(?;) - 180(402 cos(2?;) + 154 cos(3?;) - 804 cos(4?;) 
+ 973 cos(5?;) - 115 cos(6?;) - 1308 cos(7?;) + 374 cos(8?;) + 251 cos(9?;) 
+ cos(10?;) + 326 cos(ll?;) - 60 cos(12?;) - 125 cos(13?;) ) + ?;(1200(159 
cos(?;) - 968 cos(2?;) + 1401 cos(3?;) - 459 cos(4?;) + 995 cos(5?;) + 1205 
cos(6?;) - 223 cos(7?;) + 73 cos(8?;) - 44(5 cos(9i;) + 6) - 247 cos(10?;) ) 
sm^{v) + 4?;(70204?;2 + 98445) cos(2?;) + ?;((366842?;2 - 234180) cos(3?;) 
+ 30(9964 cos(4?;) + 32995 cos(5?;) - 25405 cos(6?;) - 22548 cos(7?;) + 
4066 cos(8?;) + 105 cos(9?;) + 4791 cos(10?;) + 4542 cos(ll?;) - 1516 
cos(12?;) - 1219 cos(13?;) ) + ?;(2?;(15728640?;(294 cos(?;) + 70 cos(3?;) - 
1) sin^d) cos^^f) - 517606 cos(4?;) - 235773 cos(5?;) + 227655 cos(6?;) 
+ 127788 cos(7?;) + 14656 cos(8?;) - 3946 cos(9?;) - 39071 cos(10?;) - 
18661 cos(ll?;) + 8590 cos(12?;) + 4745 cos(13?;) ) - 3(222705 sin(?;) - 
265046 sin(2?;) + 119729 sin(3?;) - 332466 sin(4i;) - 319876 sin(5?;) + 
278964 sin(6?;) + 170124 sin(7?;) - 3966 sin(8?;) + 739 sin(9?;) - 48110 
sm(10?;) - 30601 sin(ll?;) + 12050 sin(12?;) + 7780 sm(13?;) )))))) 
blnum = 17394617548800?;^ cos^^d) sin^^d) (?;(8(95239?;'' + 132?;^ - 
6660) cos(?;) + 720(37 cos(2?;) + 36 cos(3?;) - 129 cos(4?;) + 92 cos(5?;) 
+ 40 cos(6?;) - 36 cos(7?;) + 110 cos(8?;) - 30 cos(9?;) - 77 cos(10?;) + 
8 cos(ll?;) - 3 cos(12?;) + 4 cos(13?;) + 8 cos(14?;) )+?;(- 160(3540 
cos(?;) - 4081 cos(2?;) + 4258 cos(3i;) - 6420 cos(4?;) - 2416 cos(5?;) 

- 1709 cos(6?;) - 3480 cos(7?;) + 1352 cos(8?;) + 946 cos(9?;) + 420 
cos(10?;) + 452 cos(ll?;) - 122) sm^{v) + 22?;(9223?;2 - 12168) cos(2?;) 
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+ w(4(76560 - 94681^2) cosiSv) + 48(7741 cos(4w) - 4350 cos(5v) + 
1490 cos(6t;) - 4837 cos(7w) - 5698 cos(8f) + 2835 cos(9i;) + 2127 
cos(10t;) + 178 cos(ll?;) + 205 cos(12i;) - 228 cos(13?;) - 176 cos{Uv) 
) + u(2i;(22020096i;(70 cos(2?;) + 10 cos{Av) + 63) sm^{^) cos^^(|) - 
55381 cos(4u) - 74360 cos{5v) - 19668 cos(6i;) + 94024 cos(7v) + 43706 
cos(8i;) - 24996 cos(9t;) - 14833 cos(10i;) - 3270 cos(llt;) - 1099 cos(12?;) 
+ 1808 cos(13i;) + 984 cos(14t>) ) + 210232 sm{v) - 456433 sm{2v) + 
711254 sin(3t;) + 322575 sin(4t;) - 9720 sm{5v) + 108420 sm{6v) - 
311124 sm{7v) - 200610 sm{8v) + 107672 sin(9v) + 68535 sin(10t;) + 
11650 sin(llw) + 6143 sm(12t;) - 8148 sm(13u) - 5090 sm{Uv) )))) - 
4(18341?;5 + 1344^3 - 2520t; - 3360 sm{v) + 210 sm{2v) + 2100 sin(3u) 

- 3150 sin(4i;) + 2520 sin(5u) + 2520 sin(6f) - 3150 sin(7i;) + 2520 
sm{8v) + 210 sm{9v) - 3150 sin(10i;) + 420 sin(ll?;) + 210 sm(12?;) + 
420 sm(14?;) )) 

btnum = - 1 73946 17548800t;6 cos^°(|) sin^^d) (6(165878^5 - 9859^7^ - 
9930u - 630 sin(i;) + 2520 sin(2i;) + 2310 sm{Sv) - 5670 sm{4v) + 1680 
sm(5f) + 1680 sm{6v) - 5880 sin(77;) + 2310 sin(8t;) + 1680 sm{9v) - 
840 sin(10t;) + 2310 sin(llt;) - 210 sin(12v) - 840 sin(13?;) - 210 sin(15t;) 
) + v(6(5893v'^ - 52145^2 + 9930) cos(v) + 4(99685^"^ + 104103^^ - 
13590) cos(2i;) + 180( - 361 cos(3t;) + 787 cos{4v) - 740 cos(5t;) + 212 
cos(6?;) + 1104 cos(7t;) - 377 cos(8f ) - 40 cos(9t;) - 66 cos(10f ) - 395 
cos(llt;) + 65 cos(12t;) + 78 cos(13t;) + 12 cos(14i;) + 23 cos(15u) ) 
+ v(v(2(96427u2 + 513) cos(3v) + (203778 - 842996^^) cos(4u) - 6( 

- 131620 cos{5v) + 116572 cos(6f) + 82296 cos(7t;) + 145 cos(8?;) + 
16904 cos(9?;) - 28534 cos(10t;) - 23389 cos(ll?;) + 4087 cos(12u) + 
2910 cos(13't;) + 1236 cos(14i;) + 925 cos(15w) ) + t;(2i;(44040192'u(140 
cos^iv) + 3 cos(5i;) ) sin^d) cos"(§) - 147678 cos(5t;) + 123852 cos(6t;) 
+ 76254 cos(7?;) + 47602 cos(8u) + 14264 cos(9z;) - 35354 cos(10i;) - 
18078 cos(ll?;) + 3078 cos(12u) + 2109 cos(13u) + 1116 cos(14?;) + 
603 cos(15v) ) - 3(227030 sm{v) - 365012 sin(2?;) + 74524 sm{3v) - 
328361 sm{4v) - 248102 sin(5f) + 215592 sin(67;) + 107046 sm(7t;) 
+ 40423 sm{8v) + 25032 sm{9v) - 49128 sin(10i;) - 29993 sin(llt;) + 
5145 sin(12?;) + 3537 sin(13?;) + 1748 sm(14i;) + 1072 sin(15f ) ))) - 
240( - 2688 cos{v) + 6303 cos(2u) - 8687 cos(3w) + 1939 cos(4?;) - 6433 
cos(57;) - 5887 cos(6?;) + 371 cos(7t;) - 1323 cos(8f) + 1397 cos(9i;) + 
1271 cos(lOu) + 200 cos(lli;) + 207 cos(12^;) + 2440) sm^{v) ))) 
blnum = 3478923509760u^ cos^°(f ) sin^^^i^ (?;(4(2852497?;^ + 170640?;^ 

- 207720) cos{v) + (2064158f - 2922480t'2 + 223920) cos(2t;) + 720(586 
cos(3t') - 1710 cos(4t;) + 1310 cos(5w) + 560 cos(67;) - 408 cos(7t;) + 
1504 cos(8i;) - 484 cos(9u) - 912 cos(lOu) + 70 cos(llu) - 178 cos(12u) + 
74 cos(13?;) + 121 cos(14v) + 6 cos(15i;) + 11 cos(16v) ) + ?;(i;(8(485010 

- 519679t;2) cos(3?;) + 240(20686 cos(4t;) - 12039 cos(57;) + 4270 cos(6u) 

- 14162 cos(7i;) - 14606 cos(8?;) + 6470 cos(9t;) + 4182 cos(10t;) + 1427 
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cos(llt;) - 7( - 178 cos(12t;) + 93 cos(13t;) + 69 cos(14t>) + 8 cos(15f ) 
) - 41 cos(16t;) ) + t;(2t;(110100480f(210 cos(2i;) + 42 cos{4v) + 2 
cos(6f) + 175) sin^d) cos^^{^) - 585830 cos{4v) - 1144400 cos(5t;) - 
249340 cos{6v) + 1098208 cos(7?;) + 519396 cos(8i;) - 196056 cos(9t;) - 
133908 cos(10i;) - 82040 cos(llv) - 35842 cos(12w) + 22516 cos(13t;) + 
12929 cos(14?;) + 1894 cos(15f ) + 1019 cos(16f ) ) + 3(1368670 sm(?;) - 
1910495 sm{2v) + 3490840 sin(3f ) + 1317326 sin(4t;) + 149620 sin(5i;) 
+ 522816 sin(6?;) - 1381824 sm(7t;) - 805504 sm{8v) + 337776 sin(9i;) 
+ 212736 sin(10?;) + 109796 sin(lli;) + 61986 sm(12?;) - 35824 sin(13z;) 

- 22679 sin(14t;) - 3050 sm(15?;) - 1849 sm(16t;) ))) - 2400(4366 cos(t;) 

- 4746 cos(2t;) + 4156 cos(3^;) - 6806 cos(4u) - 2270 cos{5v) - 2207 
cos{6v) - 3446 cos(7i;) + 1146 cos(8?;) + 572 cos(9i;) + 557 cos(10i;) + 
544 cos(llt;) + 38 cos(12z;) + 38 cos(13t;) - 654) sm''^(t;) )) - 2(227863?;^ 
+ 369240?/^ - 1054801; - 75600 sm{v) - 6300 sm(2t;) + 56700 sin(3t;) - 
80640 sin(4i;) + 69300 sin(5u) + 70560 sin(6w) - 80640 sin(7t;) + 70560 
sin(8i;) - 5040 sin(9t;) - 80640 sin(lO'y) + 13860 sm(ll?;) - 5040 sin(12t;) 
+ 1260 sin(13f) + 13860 sm(14t>) + 1260 sm(16i;) )) 

66 = - 3478923509760i;6 cosi°(f ) sm^'^{^) (i;((468896w'^ - 2845950?;2 
+'504540) cos(u) + (3794666^"^ + 3508410?;^ - 424260) cos(2z;) + 180( 

- 3359 cos(3t;) + 6514 cos(4?;) - 5660 cos(5u) + 1975 cos{6v) + 8736 
cos(7z;) - 3088 cos(8t;) + 88 cos(9t;) - 711 cos(10t;) - 3340 cos(ll?;) + 526 
cos(12i;) + 487 cos(13?;) + 133 cos(14i;) + 238 cos(15v) + 4 cos(16t;) + 
7 cos(17?;) ) + ?;(?;(6(218182?;2 + 75025) cos(3?;) - 30( - 48938 cos(4?;) 

- 205188 cos(5t;) + 189209 cos(67;) + 123872 cos(77;) + 9360 cos(8?;) + 
34024 cos(9z;) - 47753 cos(10i;) - 37460 cos(ll?;) + 4522 cos(127;) + 2997 
cos(13t;) + 2499 cos(14i;) + 1846 cos(15t;) + 68 cos(16w) + 49 cos(17?j) 
) + t;(2^;(31457280'i;(49(25 cos{v) + 9 cos(3i;) + cos{5v) ) + cos(7i;) 
jsin^d) cos"(|) - 3244536 cos(4?;) - 1062160 cos(5i;) + 806275 cos(6w) 
+ 541016 cos(7w) + 436792 cos(8t;) + 145188 cos(9w) - 269891 cos(10z;) 

- 141460 cos(llt;) + 11496 cos(12z;) + 10442 cos(13f ) + 10763 cos(14z;) 
+ 5913 cos(15t;) + 274 cos(16t;) + 147 cos(17^;) ) - 5461475 sm{v) + 
9728100 sm(2t;) - 1651927 sm{3v) + 8272210 sin(4t;) + 5719296 sm(5t;) 

- 4891402 sm{6v) - 2271864 sm{7v) - 1253544 sm{8v) - 729384 sin(9^;) 
+ 1169646 sin(10i;) + 708816 sin(llv) - 69694 sin(12?;) - 53079 sin(13w) 

- 51572 sm(14t;) - 31759 sm(15t;) - 1350 sm(16i;) - 812 sin(17i;) )) - 400( 

- 17068 cos{v) + 33671 cos(2t;) - 46168 cos{3v) + 9411 cos{4v) - 34507 
cos(5f) - 29117 cos(6w) + 423 cos(7w) - 7921 cos(8i;) + 7213 cos(9w) 
+ 6135 cos(10t;) + 1367 cos(llt;) + 1363 cos(12t>) + 36 cos(13t;) + 35 
cos(14t;) + 14143) sm^{v) )) - 4( - 2116147?;^ + 59850t;-'^ + 134820v 
+ 5775 sm(v) - 26775 sm{2v) - 30345 sm{3v) + 67200 sin(4u) - 17640 
sin(5v) - 17745 sm{6v) + 70560 sm(7'u) - 30240 sin(8'i;) - 17640 sm{9v) 
+ 9135 sm(lOf) - 30240 sm(llt;) + 3360 sm(12?;) + 9135 sin(13i;) + 
105 sin(14t;) + 3360 sin(15?;) + 105 sm(17t;) )) 
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Figure 1. The stability region (s-v plane) of the classical Quinlan-Tremaine 14-step 
method and of the methods PF-D0,PF-D1, PF-D2, PF-D3, PF-D4, PF-D5 and 
PF-D6 (from left to right and from top to bottom) 
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Figure 2. The accuracy (digits) of the new methods compared to the classical one 
for the Schrodinger equation (E=989) 
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Figure 3. The accuracy (digits) of the new methods compared to the classical one 
for the Schrodinger equation (E=341) 
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Figure 4- The accuracy (digits) of the new methods compared to the classical one 
for the Schrodinger equation (E=163) 



qT114_JMC.tex; 15/11/2008; 10:35; p. 42 



